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Abstract. We present direct evidence for CO freeze-out in a circumstellar disk around the edge-on class I object 
CRBR 2422.8-3423, observed in the M band with VLT-ISAAC at a resolving power R ~ 10 000. The spectrum 
shows strong solid CO absorption, with a lower limit on the column density of 2.2 x 10^* cm^^. The solid CO 
column is the highest observed so far, including high-mass protostars and background field stars. Absorption by 
foreground cloud material likely accounts for only a small fraction of the total solid CO, based on the weakness 
of solid CO absorption toward nearby sources and the absence of gaseous C^*0 J = 2 ^ 1 emission 30"south. 
Gas-phase ro- vibrational CO absorption lines are also detected with a mean temperature of 50 ±10 K. The average 
gas/solid CO ratio is ~ 1 along the line of sight. For an estimated inclination of 20° ± 5°, the solid CO absorption 
originates mostly in the cold, shielded outer part of the flaring disk, consistent with the predominance of apolar 
solid CO in the spectrum and the non-detection of solid OCN~, an indicator of thermal/ultraviolet processing of 
the ice mantle. By contrast, the warm gaseous CO likely originates closer to the star. 
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1. Introduction 

Interstellar gas and dust form the basic ingredients from 
which planetary systems are built (e.g., van Dishoeck & 
Blake |1998| , Ehrenfrcund & Charnley po5o| ) . In particular, 
the icy grains can agglomerate in the cold midplanc of cir- 
cumstellar disks to form planetcsimals such as comets. In 
the cold (T < 20 K) and dense (nn = 10^ - 10^ cm'^) re- 
gions of disks, all chemical models predict a strong freeze- 
out of molecules onto grain surfaces (e.g., Aikawa et al. 
2002| ). The low molecular abundances in disks compared 
to those in dense clouds as derived from (sub)millimeter 
lines is widely considered to be indirect evidence for freeze- 
out (Dutrey, Guilloteau & Guehn |l997t Thi et al. ^OOll) . 



Observations of gaseous and solid CO have been per- 
formed for a few transitional objects from class I to class 
II that are known to posses a disk. Boogert et al. ( ^002aD 
observed L 1489 in Taurus - a large 2000 AU rotating 
disk -, but the amount of solid CO is not exceptionally 
high (~7% of gaseous CO). This may stem from the fact 
that these systems are still far from edge-on (inclination 
^ 20° — 30°) so that the line of sight docs not intersect 



the midplanc, the largest reservoir of solid CO. Shuping 
et al. ( 2001 ) found strong CO depletion toward Elias 18 
in Taurus, but both the disk stucture and its viewing 
angle are not well constrained. More promising targets 
are pre-main-sequence stars for which near-infrared im- 
ages have revealed nebul osities separated by a dark lane 
(e.g. Padgett et al. 199£ ). The lane is interpreted as the 
cold midplane of a disk seen close to edge-on where vis- 
ible and even near- infrared light are extremely extinct. 
Among such dark-lane objects, CRBR 2422.8-3423 is a red 
{H — K~A.7) low luminosity [L = 0.36 Lq, Bontemps et 
al. 2001) object surrounded by a near edge-on disk, discov- 
ered in images with the ESO Very Large Telescope (VLT) 
at 2 /im (Brandner et al. 200C ) . Its spectral energy distri- 
bution (SED) is consistent with that of a class I object or 
an edge-on class II object with strong silicate absorption 
at 9.6 /Ltm. It is located at the edge of the p Oph cloud 
complex in the core F, ~ 30" west of the infrared source 
IRS 43 and a few arcmin south-east of Elias 29 (Motte, 
Andre & Montmerle |1998D. 



* Based on observations obtained at the European Southern 
Observatory (ESO), Paranal, Chile, within the observing pro- 
gram 164.1-0605. 



This letter reports the detection of a large quantity of 
solid CO and the presence of gaseous CO in the line of 
sight of CRBR 2422.8-3423 using the ESO- VLT (§2 and 
3). Possible contamination by foreground cloud material is 
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considered in § |j, followed by a discussion on the location 
and origin of the CO gas and dust in the disk (§5). 

2. Observations 

CRBR 2422.8-3423 was observed with the ESQ VLT- 
ANTU mounted with the Infrared Spectrometer And 
Array Camera (ISAAC) on May 6, 2002. A spectrum at 
R « 10'' was obtained in the M band using a slit-width 
of 0'.'3, which matched the excellent seeing (~ 0'.'3 at 4 
fim). The integration time was 36 minutes, resulting in a 
continuum S/N of ~ 20. The spectra of BS 6084 (BlIII) 
and BS 6378 (A2.5Va), observed immediately before and 
after CRBR 2422.8-3423, were used to remove the atmo- 
spheric features and to calibrate the spectrum in flux and 
wavelength. All the reduction steps were carried out using 
an in-house data reduction package written in IDL. 

Observations of gaseous CO (sub)millimeter lines on 
source and at an offset position 30"south were obtained 
at the James Clerk Maxwell Telescope (JCMT)Q in June 
2002. The dual-polarization receiver B3 was tuned to ob- 
serve the I^CO J = 3 ^ 2 (345.796 GHz, beam 13'.'8) 
and i^CO J = 3 ^ 2 (330.587 GHz, beam 14'.'4) hues, 
and receiver A3 to obtain C^^O J = 2 ^ 1 (219.560 GHz, 
beam 22". 2). The backend was the Digital Autocorrelator 
Spectrometer, set at a resolution of ~0.15 km s~'. The 
beam efficiencies are 7ymb=0.63 (345 GHz) and 77mb=0.69 
(230 GHz). The fluxes were calibrated against the bright 
nearby source IRAS 16293-2422. The observations were 
acquired in position switching with a throw of 1800"south 
and reduced using the SPECX software. 

3. Results and analysis 

3.1. VLT-ISAAC spectrum 

A saturated absorption feature centered at ^ 4.67 fim and 
assigned to solid CO is seen in the Af-band spectrum dis- 
played in the upper panel of Fig. |l|. No broad absorption 
at 4.62 /im, usually attributed to solid OCN^ and a sign- 
post of thermal and/or ultraviolet processing, is present, 
although the limit of r < 0.05 is not strong. Unresolved 
gas-phase CO lines are also detected. A continuum which 
takes into account absorption by silicate is fitted to the 
data. The continuum substracted spectrum is shown in 
the lower panel of Fig. ^ in transmission scale. 

The best fit to the saturated solid CO band is r « 6.2, 
corresponding to A''(COico) ~ 2.2 x 10""^^ c m~^ u sing the 
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integrated band strength of Gerakines et al. ( 1995 ). This is 
the highest solid CO column density found to date, even 
compared with ice-rich high-mass source s such as NGC 
7538 IRS9 (e.g. Chiar et al. |l99|, |1995| ). The solid CO 
profile toward CRBR 2422.8-3423 can be decomposed into 



' The James Clerk Maxwell Telescope is operated by the 
Joint Astronomy Centre in Hilo, Hawaii on behalf of the 
Particle Physics and Astronomy Research Council in the 
United Kingdom, the National Research Council of Canada 
and The Netherlands Organization for Scientific Research. 
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Fig. 1. VLT-ISAAC A/band spectrum toward CRBR 
2422.8-3423 (upper panel). The i nsert shows the VLT Ks 
image taken by Brandner et al. (2000). The upper curve 
of the lower panel shows the gas-phase model, convolved 
to the observed spectral resolution and shifted in wave- 
length to account for a source heliocentric velocity of 10 
km s^'. The middle curve displays laboratorium absorp- 
tion spectra of apolar (dotted) and polar (dashed) CO ice 
mixtures. The sum of the two models including Pfund (3 
emission at 4.65 fim is overlaid on the observed spectrum 
in the bottom curve. 



a narrow saturated component, usually ascribed to apolar 
CO in a H2 0-poor matrix, and a broad red wing at 4.685 
/im, indicative of polar CO in an H2 0-rich matrix. Using 
optical constants from the Leiden database including grain 
shape effects (Ehrenfreund et al. 1997), a search for the 
best-fitting laboratory mixtures was performed. The sat- 
uration of the line prevents any unique fit, but the apolar 
CO is best matched by a mixture CO:N2:CO2=100:50:20 
at 10 K with TV ==2.0 x 10^^ cm~^. The red-wing is fit- 
ted by HaOiCO = 4:1 at 40 K with 2.6 x lO" cm^^. 
The latter column density is comparable to that found to- 
ward L 1489, which has a luminosity of 3.7 Lq. The large 
amount of apolar CO indicates low temperatures along 
the fine of sight, since it evaporates around ^ 20 K. 
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The VLT spectrum shows the presence of narrow gas- 
phase ^^CO and ^^CO ro- vibrational hnes originating from 
levels up to J = 9 (250 K above ground). Synthetic LTE 
model spectra were fitted to the observed spectrum using 
data from the HITRAN database (Rothman et al. |l99§). 
The fit parameters are the gas temperature Tex, the col- 
umn density iVgas(CO) and the velocity broadening AV, 
assumed to be smaller than 2 km s~^ from the JCMT 
data. The limited number and the high optical depth of 
the ^^CO hnes prevent a unique fit, but the best result (see 
Fig. 1) is obtained for T^x = 40 — 60 K, which is probably 
a mean between the cold and warm components along the 
line of sight. The fit includes the ^'^CO lines. The optical 
depth of the ^^CO lines is likely underestimated so that 
the gas-phase CO column density can range from 1 to 6 x 
10^^ cm~^. Adopting a mean value of 3 x 10^^ cm~^, the 
line of sight average gas/solid CO ratio is ~ 1. Comparing 
the column of gas-phase CO with that of polar solid CO 
only, the ratio drops to ~ 0.1. The latter value is com- 
parable to that found for L 1489 (0.07, see Boogert et al. 



2002a|) b ut still higher than for EUas 18 (0.01, Shuping et 
al. ^001 ). If the lines were significantly wider than 2 km 



s as found for L 1489 (~ 20 km s ), the gaseous CO 
column density would drop by an order of magnitude. 



3.2. JCMT spectra 

The C^^O J = 2 ^ 1 emission observed with the JCMT is 
shown in the upper panel of Fig. ||. Strong ^^CO J = 3 ^ 
2 (~ 12 K peak temperature) and ^^CO J = 3^2(~7K 
peak temperature) lines are also detected, but are not dis- 
played here since their profiles are similar to that of C^*0. 
C^^O J = 2 ^ 1 is detected on source only. Compared 
with single-dish CO 3-2 and 2-1 lines from disks around 
older isolat ed pre-main-sequence stars in Taurus (e.g., Thi 
et al. ^OOlD , the hnes toward CRBR 2422.8-3423 are more 
than an order of magnitude stronger and do not show 
the double-peak profile resulting from the projection of a 
disk in Keplerian rotation. The CRBR 2422.8-3423 spec- 
trum shows three peaks, two of which have velocities sim- 
ilar to those seen toward Elias 29 and interpreted as aris- 
ing from the cloud ridge in which the source is embedded 
(A^(COgas)cioud ^ 2.9 X 10^8 cm~2, Boogert et al. |2002bD . 
The total integrated C^^O 2-1 intensity on source is 13.6 
± 4 K km s~^. The corresponding ^^CO column density 



is ( 0.5-1 ) X 10^^ cm-2 assuming Tex = 15 K (cf. Motte et 



al. |1998| for p Oph clump F) and i**O/i6O=560 (Wilson 
& Rood 1994 ). For Tex = 50 K, the column densities arc 
lowered by 30%. The non-detection of C^^^O J = 2 ^ 1 
at the offset position suggests a ^^CO column density less 
than 1 X 10^^ cm^^ for Tex = 15 K , around 20 times lower. 



4. Contamination by foreground material? 

Absorption studies toward young stellar objects located in 
the p Oph cloud complex can be dominated by foreground 
cloud(s) (Boogert ct al. p002bD . The SCUBA 850Aim map 
of p Oph by Johnstone et al. (|2000D shows that CRBR 
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Fig. 2. C^^O J = 2 ^ 1 obtained using the JCMT toward 
CRBR 2422.8-3423 on source (upper spectrum) and at an 
offset position 30" south (lower spectrum). 



2422.8-3423 is indeed at the edge of the ridge which con- 
tains both Elias 29 and the nearby class I object IRS 43. 
The lack of CO emission at the offset position indicates 
however that the bulk of the gas-phase CO is located 
within 3000 AU from CRBR 2422.8-3423. The similarity 
of the gaseous CO column densities derived from the on- 
source C^^O 2-1 emission and the VLT-ISAAC infrared 
absorption may be fortuitous since the VLT-ISAAC data 
do not probe very cold CO with small (^ 1 km s~^) line 
widths. We cannot exclude that part of the gas-phase CO 
seen in the infrared arises in a more extended envelope or 
cloud, but the fact that the CO excitation temperature is 
significantly above 10 K indicates that at least some frac- 
tion must originate close to the young star in a warm part 
of the disk. 

For solid CO, there are strong arguments that most of 
the absorption must arise in the disk. The bright nearby 
(^ 30"east) source IRS 43 was observed simultaneously 
with CRBR 2422.8-3423 with VLT-ISAAC and has a CO 
ice optical depth of only 1.87 ± 0.02, corresponding to 
Af(COico)=7 X 10^^ cm^^ (Pontoppidan et al., in prep.). 
This is at least a factor of 3 lower than toward CRBR 
2422.8-3423 even though its 850 pm fiux is a factor of 1.5 
higher. Toward Elias 29, solid CO has an optical depth of 
only 0.25, most of which is believed to be located in the 
foreground clouds. Compared to the optical depth found 
toward CRBR 2422.8-3423 (rico(CO) > 6), this amount of 
foreground material can account for only an insignificant 
fraction of the observed solid CO. Finally, if clouds happen 
to lie in front of CRBR 2422.8-3423, the moderate extinc- 
tion {Ay <10) of those clouds probably prevents them to 
harbour significant amounts of solid CO (Shuping et al. 
pOOO| ). Indeed, in a A/-band survey of more than 30 young 
stellar objects, this is the deepest CO ice band observed. 



4 



W. F. Thi et al.: Detection of abundant solid CO toward CRBR 2422.8-3423 



1.2 



n r 



CO-CO E^^ = 960K 
CO-H,0 E.'' = 1740K 




Density (cm ' 



Fig. 3. Fraction of CO molecules in the solid phase at 
different temperatures as function of density. Adsorption 
onto H2O and CO ice produce different results. Th e des- 
orption energies are from Sandfor d fc A Uamandola ( 1988 ) 
for CO-H2O and Sandford et al. ( |l988|) for CO-CO. 



5. Discussion 



A simple disk model cf. Chiang & Goldreich (1999) with 
T,=3500 K and a disk radius of 250 AU has been adopted 
to investigate the location of the solid and gaseous CO 
seen in infrared absorption. Because of its higher temper- 
ature, the gaseous CO is likely not co-located with the 
bulk of the apolar solid CO, which evaporates at ~ 20 K. 
The polar solid CO can however reside in the same re- 
gion of the disk as the CO gas at 40-60 K. Assuming that 
CO is frozen out at < 20 K (apolar CO) and < 40 K 
(polar CO), the best fit to the column densities, gas/solid 
CO and polar ice/apolar ice ratio is obtained for j=20±5°. 
This inclination is consistent with the flux asymmetry seen 
in the near-infrared image of Brandner et al. ( ^000|) . For 
such line of sight, most of the CO ice is located above the 
midplane in the outer disk, whereas the CO gas is found 
in the warm inner disk. Thus, the overall CO depletion 
could be significantly higher than the ratio of ^ 1 found 
here. 

Several time-dependent chemical models were run to 
quantify the gas/solid CO ratios in different density 
and temperature regimes. The models simulate gas-phase 
chemistry, freeze-out onto grain surfaces and thermal as 
well as non-thermal evaporation. Cosmic-ray induced des- 



orption is modeled as in Hasegawa & Herbst (1993), which 
may be an overestimate for large grains in disks (Shen et 
al., in prep.). The sticking coefficient was set at 0.3 to ac- 
count for other non-thermal mechanisms (e.g., Schutte & 
Greenberg 1997) and photodesorption is assumed ineffec- 



tive. At T > Tcvap thermal desorption dominates, whereas 
at T < Tcvap cosmic-ray desorption prevails. In the model, 
Tevap is 20 K for apolar CO ice {Edcs^960 K) and 40 K 
(£;des=1740 K) for polar CO. To illustrate the eflFects of 
thermal and cosmic-ray induced desorptions. Fig. |^ shows 
the gas/solid CO ratio at chemical equilibrium. In cold 



(T < 20 K) but moderately dense (lO'' - 10^ cm^^) j.g_ 
gions of the disk around CRBR 2422.8-3423, only a smaU 
amount of CO is depleted onto grains in the apolar form, 
but much larger fractions > 50% can occur at higher den- 
sities (- 10*^ - 10^ cm-3). 

In summary, we detected a large amount of solid CO in 
the fine of sight toward CRBR 2422.8-3423. The majority 
of this ice is likely located in the flaring outer regions of the 
edge-on circumstellar disk. Very high resolution {R > 10^) 
near-infrared spectroscopy is needed to reveal the gaseous 
CO line proflles and thus their origin and the gas dynam- 
ics in the inner disk. Future submillimeter interferometer 
data can probe the velocity pattern and excitation con- 
ditions of the gas as functions of disk radius, whereas 
mid-infrared spectroscopy with, e.g., the Space Infrared 
Telescope Facility (SIRTF) will allow searches for other 
ice components, in particular solid CO2. 
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